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The reaction dynamics of axisymmetric deformed ^*Mg + '^'^Mg collisions have been investigated 
systematically by an isospin-dependent quantum molecular dynamics (IDQMD) model. It is found 
that different deformations and orientations result in apparently different properties of reaction 
dynamics. We revealed that some observables such as nuclear stopping power (R), multiplicity of 
fragments, and elliptic flow are very sensitive to the initial deformations and orientations. There 
exists an eccentricity scaling of elliptic flow in central body-body collisions with different deforma- 
tions. In addition, the tip-tip and body-body configurations turn out to be two extreme cases in 
central reaction dynamical process. 
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Aligned experiments investigating how deformed ^^^Ho 
target affects the total neutron reaction cross section 



from 2 to 125 MeV 
l5<Ea< 23 MeV 



and scattering of a particles with 
were carried out about forty years 
ago. The similar case occurs in nanoscale physics that the 
initial shape of hot droplets also has significant effects on 
fragmentation process in the molecular dynamics (MD) 
framework [3J. It is expected that deformed nuclei in- 
duced heavy-ion collisions (HICs) can result in obviously 
different properties of dynamical processes and final state 
observables compared with spherical cases. There are 
some reports about deformed U -|- U collisions at rela- 
tivistic and ultrarelativistic energies and it is suggested 
that deformed U + U collisions are more likely to cre- 
ate quark-gluon plasma (QGP) and may resolve many 
outstanding problems |4[-[lll| . The deformation effects on 
reaction cross section [1^, elliptic flow [l3| and heavy- 



ion fusion [ij, [15| was also discussed recently. On the 
other hand, polarized target and beam have been widely 
applied related with spin effects in HICs [16| especially 
for the total and differential reaction cross section mea- 



surement of aligned deformed beams such as ^Li 17 [ and 
23Na [13. 

The spin polarized beams have been greatly promoted 
by projectile- fragmentation reactions recently [l9|, which 
brings large angular momentum into fragment spin. Not 
only the fragmentation process itself produces spin po- 
larized fragments but also the produced spin orientated 
beam of deformed nuclei can provide valuable informa- 
tion on shape effects during collisions [20J. Therefore, 
it is very necessary to consider the degree of freedom of 
initial deformation since so many radioactive nuclei far 
from /3-stability line may have large deformation. How- 



ever, the knowledge about collisions induced by deformed 
nuclei is very poor especially at intermediate energy. 

Due to the distinct differences in overlap region of de- 
formed nuclei collisions, collisions of aligned deformed 
nuclei may give a clearer and deeper insight into the re- 
action mechanism such as the process of multifragmenta- 
tion and the development of collective flow. The different 
orientational collisions also have the advantage in fixing 
the uncertain behavior of density dependent symmetry 
energy, which is an elementary open problem related not 
only to many problems in nuclear physics but also to a 
number of important issues in nuclear astrophysics [2l| . 
Besides the advantage in studying reaction mechanism 
and dynamics, highly deformed nuclei induced reactions 



22| 



may also inspire exotic nuclei research such as halo 
and cluster phenomena |23|. 

In this paper, ^"'Mg -I- ^"'Mg collision system is used to 
investigate the initial deformation and orientation effects 
by a microscopic transport model: the IDQMD model 
[24|, which was developed from the quantum molecular 
dynamics (QMD) model 25[. The main advantage of 



the QMD model is that it can explicitly treat the many 
body state of collision system. So it contains correlation 
effects to all orders and can treat the fragmentation and 
fluctuation of HICs well. 

In this calculations, soft and hard nuclear equation of 
state (EOS) with the incompressibility oi K = 200 and 
380 MeV, respectively, are used for comparison. Here 



the strength of symmetry potential C. 
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and experimental parameterized nucleon-nucleon cross 
section which is energy and isospin dependent are used. 
^^Mg is approximately treated as a sharp-cutoff ellip- 
soid with large quadrupole deformation parameter: /32 
= 0.416 [26|. For comparison, systematical calculations 
for tip-tip (body-body) collisions of ^^Mg -I- ^^Mg with 
^2 = 0, 0.05, 0.1, 0.2 (all four cases with the same root 
mean-square radius) at different energies and impact pa- 
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FIG. 1: Schematic representation of tip-tip and body-body colli- 
sions. In the coordinate system, z-direction is defined as the inci- 
dent direction and the impact parameter (6) is labeled as x-axis. 
Only b to ^^Mg's long axis in body-body collisions is considered. 
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FIG. 2: (Color online) (a): R and (b): total binary collision num- 
ber (A^i,) as a function of incident energy at freeze-out time, (c): 
The ratio of non-spherical Nf, to spherical one. 

rameters are carried out . The schematic plot of tip-tip 
and body-body collisions is illustrated by Fig. [T] 

Firstly, we discuss the nuclear stopping power (R = 



-Tf 



P, 



i± 



/E- 



p. 



^11 



where A refers to the 



sum of projectile mass number and target mass number, 
Pi± = (Pi + Plyf^'^. Pi\\ = Piz in the cm. reference 
system [27|) of different orientational collisions. R can 
be used to describe the momentum dissipation and the 
degree of thermalization. Fig. [2] (a) shows that central 
body-body collisions lead to larger R than central tip- 
tip collisions below 50MeV/nucleon while the situation 
reverses when incident energies exceed 75MeV/nucleon. 
The more prolate ^"^Mg is, more obvious differences ap- 
pear. The spherical case lies between tip-tip and body- 
body collisions at all calculated energies. The larger R 
of tip-tip collisions at higher energy is in agreement with 
the resuh at 0.52GeV/nucleon by ART model 'lo'l. How- 
ever, the inversion of R between tip-tip and body-body 
collisions is first observed. It reflects the different roles 
of the initial space configurations vs. energies. 

It is known that the reaction mechanism at intermedi- 
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FIG. 3: (Color online) (a), (b) and (c) represent the time evolu- 
tions of R, maximal density Dmax and Ni, in central collisions(b = 
Ofm), respectively. The long-dashed, short-dashed and dash- 
dotted lines are drawn to mark the characteristic time of the colli- 
sions. The time structure of the Ni, is synchronous with the den- 
sity evolution. Time evolution of R show tip-tip, body-body and 
sphere-sphere collisions experience different touching, compressing 
and expanding processes from t = to freeze-out stage. 



ate energy is dominated by mean field, binary collisions 
and Pauli blocking. Since the IDQMD model can treat 
the three components explicitly, it is very convenient to 
find out the factors which dominate the stopping power 
at different energies. As represented in Fig. [2] (b) and 
(c), tip-tip collision numbers are higher than body-body 
ones at all considered energies. It means that binary 
collision cannot be responsible for the inversion of R, 
while the mean field must play a very important role. 
Fig. [3] shows how mean field and binary collision take 
effect in dynamical process. The peak of density corre- 
sponds to the most intensive stopping process but the 
R has not reached maximum. The departure between 
projectile-like and target-like continues contributing the 
nuclear stopping power. Through the different stopping 
behaviors of tip-tip and body-body collisions vs. ener- 
gies, the time evolutions of R show that when reaction 
proceeds more quickly, the larger stopping power can be 
achieved. So the stopping power can be regarded as a 
measurement of time scale of dynamical process as well 
as an observable of momentum dissipation. 

Due to the larger projectile-target overlap region, 
body-body configurations build up stronger mean field, 
which lead to more violent one-body scattering. How- 
ever, transparency effect of the nuclear medium becomes 
more and more important when incident energy rises. 
The tip-tip configurations are less transparent, which 
leads to larger two-body collisions and stopping power. 
Therefore, the underlying mechanics of the inversion 
of R between tip-tip and body-body collisions is that 
body-body configurations build up stronger mean field 
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FIG. 4: (Color online) Energy dependence of fragment multiplic- 
ity of tip-tip, body-body and sphere-sphere collisions at different 
reduced impact parameters (ft^ed = b/bmax, where bmax refers to 
the maximal impact parameter for different cases). 

at lower energies, where one-body scattering is predomi- 
nant. Whereas two-body coUisions become more impor- 
tant in tip-tip configurations at higher energies. 

Since the IDQMD model can treat fragmentation of 



hot nuclei [28|, |29[ well. It is appropriate to investigate 
the fragmentation observables. As shown in Fig. 21 the 
fragment multiplicity has strong correlation with stop- 
ping power. Body-body collisions have minimal multi- 
plicity at all impact parameters at higher energies while 
tip-tip collisions have the maximal one. So this behavior 
is consistent with that of stopping power at higher ener- 
gies. It can also be seen from charge distributions in Fig. 
[5] that the tip-tip and body-body collisions are two ex- 
treme cases and the sphere-sphere collisions lie between 
them. Therefore, the fragment observables also confirm 
the similar picture indicated by R. 

The Body-body collisions with b — fm will produce 
large collective motions due to the different initial ge- 
ometry from spherical nuclei. Anisotropic flow method 
has been developed to measure the anisotropy of parti- 
cle momentum space which related to the EOS and nu- 
clear reaction dynamics 301432 1 . The azimuthal distribu- 
tion of fragments can be expressed by Fourier expansion 
[33| ^ oc 1 -f 2^^^ v„cos(n(/()), where is azimuthal 
angle between the transverse momentum of the particle 
and the reaction plane. The coefficient f„ is defined as 
anisotropic flow parameter, among which V2 denotes el- 
liptic flow. It can be calculated in terms of single-particle 

2 _ 2 

averages: V2 = (cos (20)) — { ^'Ij^l ). Nucleon's V2 in- 
duced by deformed U -I- U collisions has been studied 
by ART model [5|, llO| and optical Glauber model [13| at 
relativistic energies recently. It seems that the most cen- 
tral body-body collisions give rise to largest V2 because of 
the strongest shadowing effect in the reaction plane [5|. 
Thus V2 of central body-body collisions are most appro- 
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FIG. 5: (Color online) Charge distributions of tip-tip, body-body 
and sphere-sphere collisions at different b,,^^ and incident energies. 



priate for investigating the EOS. However, V2 developed 
from deformed nuclei collisions is unknown at intermedi- 
ate energy and it's interesting to study their deformation 
and orientation effects. 

V2 of light fragments are shown in Fig. [SJ in which the 
eccentricity (e) is calculated by maximal geometry over- 
lap region: e = Y^i^l - vDl Y.'kA + vT)- Cental tip-tip 

i i 

and sphere-sphere collisions do not have obvious V2 be- 
cause of the transverse symmetry of overlap region while 
the V2 of cental body-body collisions has non-zero value. 
The negative sign of V2 at higher energies is in agreement 
with deformed U -I- U collisions by ART model [5|, [lO| . 
The positive V2 at lower energies and the alteration of 
sign for V2 are first observed in central body-body col- 
lisions. At higher energies the violent two-body colli- 
sions in overlap region build the anisotropy pressure and 
it prompts fragments emission from in-plane preferential 
to out-of-plane preferential. The heavier fragments have 
larger V2 , which is consistent with ref. [3J] . V2 of central 
body-body coUisions (/Sa = 0.05, 0.1, 0.2, 0.416) can be 
scaled together by e from low energies to high energies. 
While scaled by the same e amplitude as the deformed 
^''Mg collisions, V2 for mid-central spherical ^^Mg colli- 
sions shows different behaviors especially for higher en- 
ergies. Therefore, the scaling of V2 indicates that the 
geometric shapes of participants play an essential role in 
collective flow of central body-body collisions. 

The energy excitation function of V2 at mid- 
central sphere-sphere collisions varies from positive (in- 
plane, rotational-like emission) to negative (out-of-plane, 
"squeeze-out" pattern) [3l|, |35| . This energy point is so- 
called transition energy, which is near lOOMeV/nucleon 
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FIG. 6: (Color online) V2 excitation function of light fragments 
at mid-rapidity (—0.5 < Y < 0.5) of deformed and spherical col- 
lisions. The dashed lines are drawn to guide the eye. Left: V2 in 
central collisions with fe = fm; Right: scaled V2 with eccentricity 
e in central body-body collisions and non-central spherical ^''Mg 
collisions. The spherical ^''Mg collisions with b = 1.45 fm and 
2.85 fm have the same absolute value of e as the deformed central 
^''Mg collisions with 1^2 = 0.2 and 0.416 , respectively. 
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FIG. 7: (Color online) Average V2 excitation function of light 
fragments at mid-rapidity (—0.5 < Y < 0.5) for deformed and 
spherical collisions with soft and hard EOS. V2 is averaged with b 
from to bmax for body-body and sphere-sphere collisions. 



3l| . For spherical collision system, there exist three 
competing components affecting the transition energy: 
(1) rotation of the compound system, (2) expansion of 
the hot and compressed participant zone, (3) shadow- 
ing of the colder spectator region [31[ . Only the expan- 
sion survives in central spherical collisionsfSB', 'ST"], which 
merely generates azimuthal symmetric flow. However, 
central body-body collisions have bulk transverse asym- 
metry overlap region and there is no rotation effect. Also 
the shadowing is different from mid-central collisions of 
spherical nuclei. Therefore, it provides an ideal tool to 
understand how the azimuthal pressure, expansion and 
flow development from the almond-shape overlap, which 
are all related with the extraction of the EOS. Average 
V2 is shown in Fig. [7] with soft and hard EOS. Hard 
EOS enhances V2 for both spherical and deformed colli- 
sions. Deformed configuration gives rise to larger V2 than 
sphere-sphere configuration for both soft and hard EOS. 
In summary, deformed ^'*Mg -I- ^^Mg collisions have 
been studied systematically by IDQMD model. The in- 
version of R vs. energies between tip-tip and body-body 



collisions reflects the two different conflgurations play dif- 
ferent roles on reaction dynamics. The fragment observ- 
ables also show different behaviors for the two configu- 
rations. The sphere-sphere collisions lie between tip-tip 
and body-body collisions in nuclear stopping and frag- 
mentation. Moreover, the excitation functions of V2 for 
different deformed central body-body collisions can be 
scaled on a similar curve by eccentricity. V2 averaged by 
impact parameter (collision configuration is represented 
by Fig. [1]) in deformed collisions is stronger than that of 
spherical collisions for both soft and hard EOS. The large 
V2 developed from cental body-body collisions have ad- 
vantages in studying the EOS and transition energy. Tip- 
tip collisions can be used to study the liquid-gas phase 
transition in finite nuclear systems due to the longer col- 
lision time. In addition, deformed nuclei collisions may 
have some implications on halo and cluster structure re- 
search. Therefore, the merits of collisions of deformed 
nuclei can shed light on the studies of both the nuclear 
structure and the reaction dynamics from low energies to 
relativistic energies. 
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